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Juan Cal VO* Matteo Novaga^ Giandomenico Orlandi^ 

Dedicated to the memory of Vicent Caselles 


Abstract 

We show existence and uniqueness results for nonlinear parabolic equations in noncylindrical 
domains with possible jumps in the time variable. 
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1 Introduction 

In recent years there has been a renewed interest in problems related with partial differential equations 
formulated in domains that change over time. This is partly due to the fact that a number of problems 
in mathematical biology are naturally posed on growing domains (e.g. developing organisms or prolifer¬ 
ating cells, see for instance [I3l|22l[20]) or domains that evolve in some particular way. Such issues have 
originated a wide amount of mathematical research, let us mention [291 El [a US]. To this we should 
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add more classical engineering applications like fluids or gases in settings as channels or pipes with 
confining walls that may be displaced, removed or brought in at will. A sample of different applications 
of partial differential equations in evolving domains can be found in the recent survey paper [2T]. In 
fact, it is not hard to conjecture that new applications will involve equations in moving domains in the 
future. Apart from that, partial differential equations posed on non-cylindrical domains are interesting 
also from the purely mathematical point of view. 

This has led to an outburst of works on this subject in the literature that added up to some classical 
works [a ESI [ini [Ml El IE] to the extent that the number of current references is overwhelming. Let us 
comment on this literature according to the approach, the assumptions on the evolution of the domain 
where the equation is posed and the types of equations considered. Many authors used semigroup 
methods to tackle these problems (see for instance mm and references therein), but other approaches 
include adding a time viscosity [9], mapping the spacetime domain to a cylindrical domain [3] or use 
De Giorgi’s minimizing movements mm- As regards time variations of the domain, it is customary to 
impose some sort of continuity (for instance Lipschitz continuity [30], relaxed to Holder continuity in [9] 
and to absolute continuity in [2l|), alternatively a monotonicity condition can be used (i.e. expanding 
domains diE]) or Reinfenberg-type domains can be considered m- Concerning the type of equations, 
most of the works focus on parabolic equations which are assumed to be linear or in divergence form 
(see however [SIEZIISIEH] where also other operators are admitted). 

In this paper we are interested in well-posedness of parabolic equations in divergence form, in 
bounded domains that evolve in time. More precisely, we deal with the Cauchy-Dirichlet problem, in 
a formulation that allows boundary conditions to depend on time. 

Let us discuss what are the novelties of this work with respect to the already existing literature. 
First, we introduce a simple approach to construct solutions, which consists in performing a time 
slicing of the domain, and then solve a family of approximating equations in cylindrical domains. The 
simplicity of this approach may allow to use it as a starting point for devising numerical methods for 
this sort of problems. Despite its simplicity, we are not aware of other works where such a slicing 
strategy is used. Our approach allows to deal with nonlinear equations, which include the parabolic 
p-Laplacian as a particular case. Also, our slicing technique applies to quite general variations on 
the domain over time: we only require them to be of bounded variation, allowing for sudden jumps 
(expansions or contractions) of the domain. In particular, we do not impose any constraint on the 
topology of the evolving domains, which may differ from that of the initial domain. We are also able 
to prove uniqueness under some additional constraints on the domain (see Section [5]) . 

To our best knowledge, this generality has not been previously achieved in the literature, except for 
the case of purely expanding domains [7|. However, in [28j F. Paronetto proposes a different approach, 
which can be extended to cover quite general operators and boundary conditions. 

Possible extensions. Since our main goal is presenting a method to tackle parabolic equations in 
moving domains, we did not focus on looking for the most general possible result. For instance, for the 
sake of simplicity we chose to deal only with bounded initial data. We stress that our main idea is to 
use a time slicing to approximate the original problem by a sequence of problems defined on cylindrical 
domains. As we do not focus on any particular equation, we chose to use abstract Lions’ theory to 
provide existence for the approximating problems. However, we could also use other theories as starting 
point to provide existence of approximate solutions. If we are interested in a particular equation (the 
p-Laplace equation, say) then we will likely be using specific existence results to set up our method, 
and those will provide a much more accurate framework for the admisible set of initial conditions. 

In that line of thought, the fact that our present formulation does not allow to deal with degenerate 
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equations, such as the porous media equation and its variants, could appear as a drawback. Again, we 
argue that suitable modifications of the method here proposed would allow to tackle these problems. In 
fact, even sticking to Lions’ theory, porous media equation and related ones can be treated by making 
use of the compactness results by Dubinskii carefully adapting our arguments in order to cope 
with that (see [251 Chapter I, 12]). We did not pursue this line here in order to keep the presentation 
as simple as possible. 

We also point out that we cannot deal with operators with linear growth such as the total variation 
flow or the parabolic minimal surface equation (see [8] for some results in this direction in the one¬ 
dimensional case). This is another challenging line to explore. Finally, following the same approach it 
should be possible to consider similar evolution equations on manifolds evolving in time and/or nonlocal 
operators (see mm and references therein). 
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2 Standing assumptions and main results 

Our purpose is to prove existence and uniqueness results for nonlinear parabolic equations with time 
dependent coefficients in time dependent domains. More precisely, given an open set Cl C [0,T] x 
we shall consider the following problem: 

ut(t, x) = div (A(t,x,u,Vu)) in D, 

< u(0,x)=uo(x) in D(0), (2-1) 

, u{t,x) ='tp{t,x) in9/Du9_iD, 

where we let be the outer unit normal to dCl, D(0) is the initial domain dehned in 

Assumption 12.21 and we set 

:= {{t, x) £ dQ : t > 0,1't = ±1}, 

diil := {{t,x) G dQ : {vtl < 1} = {(L^;) G dQ : |z^a;| > 0}. 

In order to establish existence and uniqueness of solutions, we shall make suitable assumptions on 
the flux vector held A, on the data uq, "0 and on the domain D. 

Assumption 2.1. The set Cl C (0, T) x is a bounded open set with Lipschitz boundary, and we let 

D(t) := {x G : {t,x) G D} t £ (0,T). 
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Note that 0(i) is an open set, possibly empty, for all t € (0,T). 


Notice that ri(t) has Lipschitz boundary for a.e. t € (0,r), and there exist the limits 

n(t±) := lim n{s) for all t G [0,r], (2-2) 

where the limit is taken in the Hausdorff topology. 

Assumption 2.2. The set n(0) := n(0+) is open and has Lipschitz boundary. 

Next we describe our assumptions on the operator A. Let Qq be an open set of such that 


Ute[o,T]^(^) CC Qq -where by CC we mean that the inclusion is compact- and let Qt '■= (0,r) x Qq. 
We shall denote by }A{Qt) the space of all Radon measures on Qt. 

Assumption 2.3. The function A : Qt x ]R x JR'^ — JR'^ is a Caratheodory map satisfying 

\A{t,x,z,f)\ < c\f,\^~^+ b{t,x), c> 0, b e Lp'{Qt), I < p < oo, ^ ^ = 1, (2.3) 

A{t,x,z,f,) ■ f > a\f\P — d(t,x), a > 0, d € L^{Qt), (2.4) 

{A{t, X, z, f) - A{t, X, z, C)) ■ if. - D > 0, (2.5) 

for a.e. {t,x) € Qt , and for all z ^ IR, ^ Moreover, we assume that 

\A{t,x,z,^) - A{s,y,w,^)\ < (a;(|t - s| + |x - y\) + C\z - u;|)|^|P“\ (2.6) 

where u is a modulus of continuity and C > 0. We assume also that 

A{t,x,z,0) =0 yz € M, a.e. {t,x) G Qt- (2-7) 

Note that (j2.5l) and (j2.7|) imply that 

A{t, X, z,f,) ■ f, > 0, a.e. in Qt, and for all z G fR, ^ G iR'^. (2-8) 

We will consider the problem (12.ip with initial and boundary conditions 

u(0,x) = uo(a^) € 7y°°(n(0)), (2-9) 

u{t,x) ='il){t,x), (t, x) G 9/0 U 9-1^, t > 0. (2-10) 

Assumption 2.4. We assume that 

G C(Qr) n LP{0, T; <’^(Qo)), (2.11) 


and 

V’t G L\Qt) n LP'{0,T-W-^’P'{Qo)). 

Let us now define the space after which we model the solutions of our problem. 
Definition 2.5. Let V be the closure of C^{Ll) with respect to the norm 

\ i/p 


( 2 . 12 ) 


|u||v := 


:= \Vv\Pdxdt^ ^ , u G C.^O). 

Notice that functions in V do not necessarily have zero trace on d±iLl or on 0(0). 
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Our concept of solution will be the following: 

Definition 2.6. We say that a function u G is a weak solution of i2.1\) if the following statements 

hold: 


1 . 

2 . 

3. 

4- 


u — Ip and A{t, x, u, Vu) G L^' (Q). 

ut G V* (note that this implies that u has a trace on and on n(0)J. 

u(0) = uq a.e. on 0(0) and u = ip a.e. on every relatively open subset of 
The following integral formulation 

rT r r rT f 


ucpt dxdt — / uo(p{0) dx + 


A{t, X, u, Vu) ■ Vpdxdt = 0 


(2.13) 


JO Jnp) Jq.{o) jo 

holds for all (p € ^([0,7") x Qq) with supp(/> CCit. 

Let us state the main existence result of this paper. 

Theorem 2.7. Let Assumvtions \2.I[]2.4\ be satisfied. Then there exists a weak solution of (12.ip in the 
sense of Definition 12.61 

Following [281 Assumption H.2], we introduce an additional assumption on the domain which we 
will need in the uniqueness proof. 

Assumption 2.8. For everyth G [0,T], there exist an open neighborhood U of to and a family of maps 
G{-,t) : Ll(to) Ll{t), with t G U Cl [0,T], such that 

- G{-,t) is a bijection for almost every t G U (1 [0,T]; 

- G{-,t) is Lipschitz continuous with its inverse for every t gU r\ [0,T]; 

- G{x,-) and |VG(x,-)| are absolutely continuous for almost every x G O(to); 


- |VG(-,t)| G L^(0(to)) for every t G U r\[0,T] and dt\VG\ G L^(0(to) x (0,T)). 


Note that this assumption does not allow ’’jumps” of the sections Lift). However, we could work 
in a more general framework in which the conditions in Assumption 12.81 break down for a finite set of 
times; we comment on this in Remark 15.31 below. 

Let us state now our uniqueness result. 

Theorem 2.9. Let Assumptions \2.1\\2.4\ and \2.8\ be satisfied. Then the solution of (12. ip is unique in 
the class of weak solutions. 


3 Construction of approximate solutions 

Let us divide the interval [0,T] into sub-intervals 0 = to < ti < ... < t^-i < t^ = T. The points 
are chosen so that: 

1. D{ti) has Lipschitz boundary for all i G {0,..., A — 1}, 

2. (|2.3l) ~ (l2.7p hold for a.e. x G Ll{ti) and for all z G M, f G 


5 




3. ti are Lebesgue points of V’(i) G LP{0,T; Wq'^{Qq)) and S W^'^{Qq), 

4. ti are Lebesgue points of the map t € [0, T] —>■ A{t) G L^(Qo x {—R, R) x -B(0, R)) for any R > 0, 
being B(0, R) the open ball centered at zero with radius R, 

5. A := maxfc=o,...,Af-i \tk — tk+i \ —)• 0 as iV —>■ oo. 

Let Ik = [tk,tk+i)- We iteratively solve the parabolic problem 

Uj = div ^A{tk,x,u^,Vu^)^ , t ^ Iki X ^ 

u^{t,x) = 'ijj{t,x), t ^ Ik, X G dll{tk) 

(3.1) 


x), X G fl(tk) n fl(tk-i) 

-ip(tk,x}, X e f}(tk)\fi(tk-i). 

If to = 0 we let tt'^(0,x) = uo(x). Notice that the iterative initial condition for t = tk makes sense 
thanks to the continuity properties of u^~^, see (133]). 


3.1 Study of the model problem on a time slice 

Let flo be an open bounded set in with Lipschitz boundary. Let A{x,z,^) be such that ()2.3p - ()2.7p 
hold a.e. in x G flo and for all 2 G IR, ^ G Let us consider the problem 


ttt = div (d(x, u, Vtt)) t G [0, T], X G rio) 

< u{t,x) ='ip{t,x) t G [0,T], X G drio) (3-2) 

u(0, x) = uo(x) X G no¬ 


where Ip satishes (I2.1ip - (|2.12p and no G L‘^{Qo). 

Definition 3.1. We say that a function u G L^((0,T) x Qq) is a weak solution of if u € 

LP{0,T-,W^^P{no)), d(x,n,Vn) G Lp'{{0,T) x Ho)), 



upt dxdt 



uo4>{0) dx + 



A{x, u, Vn) • Vpdxdt = 0 


holds for all p G I1([0,T) x Ho), and 


(3.3) 


u{t) — V^(t) G ILQ^’^(no) a.e. t G (0,T). 

Note that, by ()2.3I) . if n G LP{{),T-,W^'P{IIq)), then A(x,u,Vu) G RP ((0,T) x Ho)). 
Proposition 3.2. Problem (|3.2D admits a unique weak solution in the sense of Definition \S. A 
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Proof. The proof is a standard application of the theory developed in [25l [26] ; we include it for com¬ 
pleteness. We consider the auxiliary problem 


vt - div (^A(t, X, V, Vv)"j = -ift 
' v{t, x) = 0 

_ v{0,x) = uo(x) - 'tp(0,x) 


t € [o,r], X G rio, 

t G [0,T], X G 9 ^ 0 , 


X G i2o- 


(3.4) 


Here 

A{t, X, z, := A(x, z -b ipit, x),^ + Viplt, x)). 

According to the notation in [25l[26], we let ff = L^(Ho), 

r if P>2, 

[ Wo^’P(Qo)nL^(Qo) if 1<P<2, 
and F = T^(0, T‘,B), so that B is dense in H and 

-divA :F^F' = Fp\o, T; B') and V’t € F', 

with 

r W-i’P'(Ho) if P>2, 

B'=\ 

[ + A2(Ho) if 1 < p < 2. 

Observe that, by our assumptions on A(x, z,f,) and if, A{t, x,z,f,) is a Leray-Lions operator (see [25l[26] l. 
Indeed, the monotonicity requirement is satisfied thanks to (12.5p . The coercivity condition follows from 
(|2.4p and Poincare’s inequality in a standard way (lower order terms are estimated thanks to (12.31) 1. 
Then thanks to Lions’ theory there exists some v G F solving (13.411 in F' . In fact, this solution verihes 
that 

vgFP{0,T,B) and vt G IP'{0,T-, B'). 

Notice that v G C(0,T; L^(Qq)) thanks to Lemma 13.31 below. 

Let now u = v + ijj. Then tt is a weak solution of (|3.2p with initial condition u{0) = uq. Clearly 
u G LP{0,T;W^’P{no)), Ut G LP'{0,T;B') and 

uGC{0,T-L'^{no)). (3.5) 

To prove uniqueness let u, v be two different solutions. Note that u — v G F. If A does not depend 
on u, we multiply the equation for {u — v)t hy u — v and integrate by parts (see e.g. (HU Chapter III]). 
Recalling (|2.5I) we have that 

= (div(A(x, Vm) - A(x, Vu)), u - 
= Vu) - A(x, Vu), Vu - 


< 


0. 




Hence ||m — f II 2 is nonincreasing and uniqueness follows. 

For the general case, consider 5 > 0 and let 

( s ii — 5 < s < 5, 

Ts{s) = if s < —6, 

[ S ifs>(5. 

Clearly Ts{u — v) € F and again after multiplication of the equation for {u — v)t by Ts{u — v)/5 and 
integration by parts we obtain that 

1 1 

- V, Ts{u - v)/6)h = --^{VTsiu - v),A{x, u, Vu) - A{x, v, Vu)) (Qo) 

< -^((Vu - Vv)x{\u-v\<s},Mx, u, Vu) - A{x, V, Vu)) lp{Qq)-lp'( no)- 
Then, using (|2.6p we get 

-■^{u — v,Ts{u — v)/6)h<^ [ |u —u||Vm —Vu||Vu|^“^dx < C [ |Vu —Vu||Vu|^“^dx. 

dt 0 {\'U‘—'v\<S} 

The term on the far right converges to zero when d ^ 0, since the integrand is in L^(Ho) and V(u—u) = 0 
a.e. where u — u = 0. As Ts/6{u — v) ^ sign(u — u), we get that ||u — u||i is non-increasing, and we 
conclude the proof. □ 

The following continuity result is standard (see for instance [25l Ch. 2, Rem. 1.2] or [3T|1. 

Lemma 3.3. Let V be a reflexive Banach space with dual V. Let H be a Hilbert space that we identify 
with its dual. Assume that V C H C V' with the injection V d H being dense. Then, u G LP(0,T;1/) 
together with ut € IT {0,T-,V') imply that there is a representative of u which is continuous from [0,T] 
to H. 

Since € C{tkflk+i', L‘^{Ll{tk))) we can define the traces 

u^{tk+) ■= lim u^{t) u^{tk+i — ) ■= lim u^ft), 

t—t—— 

where the limit is taken in L‘^{Ll{tk)). 

3.2 The approximate solutions 

We now let 

:= {{t,x) : t G [tkflk+i], x G H(4), A; = 0,..., - 1} 

= Ufc=i,...,7V-i[4, 4 + 1 ) X H(tfc). 

Notice that does not depend only on A = maxfc=o,...,A''-i \tk — tk+i\, but depends on the entire 
sequence {4}fc- 

Lemma 3.4. converges to H in the Hausdorjf sense. As a consequence Xqa strongly in 

L^{Qt) (henee in LP{Qt) for all p < 00 ). 

Proof. The Hausdorff convergence of to H can be easily verified when H is a polyhedron. The 
claim follows by approximating a generic Ll with Lipschitz boundary with polyhedra, in the topology 
generated by the Hausdorff distance. □ 


We now glue the solutions u^{t,x) of (13.ip together and define the approximate solutions 


N-l 

u^{t,x) := (3.6) 

k=0 

N-l 

u^{t,x) := '^X[t^,tk+r)it)i^’‘it^x)xn(tk)ix)+'iPit,x)xQo\n{t^,){x)), (3.7) 

k=0 

for {t,x) G Qt- When we write u^{t,x)xn{tk)i^) ™ above formulae we intend the function which 
coincides with u^{t,x) in fl(tfc) and it is equal to zero outside fl(tfc). 

In the sequel we shall prove the compactness of and as A ^ 0. 


3.3 Estimates on 

We now derive some estimates on the approximate solutions defined in (13.6p . 

Lemma 3.5. Assume that HV'Hcxd, ||?.io||oo < C for some C > 0. Then ||u^||£,oo(qA) < C for any t > 0. 

Proof. It is enough to prove the estimate in (0,ti) x n(0). Let [•]+ denote the positive part (resp. [•]“ 
the negative part) and let C > ||V’||oo- Then the pairing of [u — C\^ with makes sense; multiplying 
(1321) by \u — C\^ and integrating by parts we get to 


2 dt 


f {[u^{t) — dx = [ — C]^divj4(0, X, Vtt^) 

Jq(o) Jn{o) 

= - [ A(0, X, u^,Vu^)V{[u^ - C]+) dx. 

Jn(o) 


dx 


There are no boundary terms present thanks to our choice of C. Note that V([rt^—C']"'') = X{u^>c}'^'^^ ^ 
so that we can use (12. 8p to ensure that the time derivative above is nonpositive. Hence, 



dx < 



([no - C\^f dx. 


Thus, if no < C* then u^{t) < C too for any t G [0, ti). This works in the same way for the time derivative 
of the integral of ([n^ + C]~)‘^ , with inequalities reversed. If we now choose C = max{||no||oo, ||V^||oo}, 
we deduce that ||n^(t)||oo < Cl. □ 


Lemma 3.6. There holds 

Pz} Pk+i r 

Y, / \Vu^it)\Pdxdt<C, 

for some constant C > 0 depending only on XI, on if and on the structural constants in Assumption 

23 


Proof. We fix k and notice that the pairing of n^ — if with n^ on (tfcj^fc+i) ^{tk) makes sense. After 
integration by parts we get 


2 dt 



— if)'^ dx 



V(n^ 


if)A{tk,x,u^,SIu^) dx — f (n^ 

Ju{tk) 


if)ift dx. 
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Notice that the last term is well-defined thanks to our assumptions on ip and to Lemma [3.5l Integrating 
the former equality on we obtain 




{u^{tk+i) -ipf dx = ^ f {u^{tk) -ipf dx- [ [ {u^ - dxi 

^ Jtk J^(tk) 

V'ijjA{tk,x,u^,Vu^) dxdt — f f Vu^A{tk,x,u^,'Vu^) dxdt 

Jo Jn(tk) 


n{tk) 

rik+i r 
tk J^dk) 


'Q{tk) 

= :/ + // + III + IV. 

Let us now control the last three terms. The second one can be easily estimated as 

r*fe+i 


II < 2C 


r^k+i r 

/ / I'lptl dxdt, (7 := max{||V^||oo, lluolloo}- 

Jtk J^dk) 


IV < — / / a\S/u^\^ dxdt + / / \d{t,x)\ dxdt. 

Jtlr Jtk JQ,(tk) 


Concerning the fourth term, using ()2.4I) we get 

rtk+i 

hk J^dk) 

In a similar way, using ()2.3I) we obtain 

C^k + l 

hk J^dk) 

Let us estimate A and B. For that we use Young’s inequality with weights: 

ePaP bP' 


III < 


/* /*£fc+i r 

/ / clVV'llVu^l^”^ dxdt-|- / / \S/p’\b(t,x) 

Jtk JfKtu) Jtk 


dxdt = A-I B. 


Then 


and 


ab < -1— - — y, e > 0,being p,p' given by (j2.3p . 

p p'eP 


B < ^\\'^'P\fLP([tk,tk+i]xn{tk)) + {[tk,tk+i]xfi{tk)) 


ce 


^-p ll^^llLP([tfe,tfe+i]xO(4)) + p/fP' \\'^^^\\LP([tk,tk+i]xnitk)) 
for any e > 0. Let us choose e so that cj{p'eP') = aj^. Collecting all the estimates, we obtain 

t; [ {u‘^{tk+i)-ipf dx+ ^ [ f \Vu^\P dxdt 

^ J^(tk) ^ Jtk JQ(tk) 


< 


/n{tk) 


'fi{tk) 


tk I t^dk) 


p 

tk+1 


1 1 /■''+! I 

+p\\^^\\U[tk,tk+i]xn{tk)) + ^WKp'iitkMxmtk)) + 


dxdt. 
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By summing up the previous inequalities from A; = 0tofe = A^ — 1, we get 


/ (u^(tAr) - ^ V / / \Vu^{t)f’ dxdt 


- ^ 


Tjf d Lr ([ifc jifc+i] Xr2(tfc)) 

fc =0 


Af-1 


'n(o) 


V 


^ p \ N — l 

2c\y^ 




ap' 


liP([ife,ife+l]xn( 4 )) 


fc =0 


+ 


1 ptu^i /• 1 /*tfc_i_i 

2 ^E/ / IV’t|dxdt+^/ 

^_Q tk " tk 


Q ^dk) 


\d{t, x)| dxdt. 


On the aid of Lemma 13.51 the thesis follows. 


□ 


Recalling the dehnition of and the assumptions on ip, from Lemma 13.61 we obtain the following 
result: 

Corollary 3.7. There exists C > Q depending only on fl, on ip and on the structural constants in 
Assum,ption \2.,‘A such that 

\\u‘^\\lp{0,T-,W^-p{Qo)) ^ ^ ■ 

In particular, the sequence {fi^} is weakly relatively compact in L^(0, T; VL^’^((5o))- 

3.4 Time compactness of 

We now show a stronger compactness property of u^. For this aim, we need the following result, proved 
in [32]. 

Theorem 3.8. Let X, B, Y be three Banach spaces such that X G B C Y. Assume that X is compactly 
embedded in B and 

F is a bounded set in L^(0,T;X), (3-8) 


\\Thf - f\\Ldo,T-h;Y) 0 as /i 0, uniformly for / € F, 
where {Thf ){t) = f{t + h) for h > 0. Then F is relatively compact in F^(0,T; B). 

Let 

N-l 

ip^{t,x) := '^X[tk,tk+A^'>^)'4’{'t^x)xQo\n{tk){x) = 'ipit,x)xQ^\n^it,x), 


(3.9) 


k=0 

-.Au _ „.Ai 


SO that we have u^{t, x) = u‘^{t, x) + ip^{t, x). 

Lemma 3.9. Let 0 < k < N be fixed. Then u^{t)xn(tk) € IT {tk,tk+i, (Llptk))) and the following 

estimate holds: 


Wt{t)Xn{tk)\\Lp' (tk,tk+i,W-^Y {Q.{tk))) ^ ^ll^^llz,P(4,4+i,lVi'P(n(tfc))) + \\^\\LP(tk,tk+i,Lp' {Q.{tk))y 
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Proof. We show the estimate by duality, 
compute 


Define Bk : 



LP(tfc,4+i, Wq’^(D(4 ))) and let (f € Bk- We 
A{tk,x,u^,'Vu^{t)) ■ Vcfdxdt. 


Hence using (12. 


f^k +1 r 

{utXn{tk),(p)B'-B^, < / (c|Vu^(t)|^’"^+ 6(t,x))|V(/>|da:dt 

Jtk d^dk) 

ftk+l / A _1 \ 

J tu 


The result follows. 


□ 


Lemma 3.10. The sequence {u^} is relatively compact in Lj^^{Q). 

Proof. We consider a cylinder C := [ti, t 2 ] x K CC D. We want to apply Theorem 13.81 with f = u^\c = 
u^\c, X = B = L^{K) and Y = W~^^'p'{K) + L^{K). Here T is a Banach space equipped 

with the norm 

||y||y := inf{||yi||^^_iy(^) + \\y2\\LpK) ■ 2/1 + 2/2 = 2/}- 

Then X C B C Y and X is compactly embedded in B. 

Notice that, since C CC H, we have 


N-l 

ut\c = Ilf Ic = X[tk,tk+i)i^)^t{t^x)xn{tk)ix)\c for N large enough. 
k=0 

Estimate (j3.8p directly follows from Lemma 13.51 In order to prove (13.9p , we notice that (with a 
slight abuse of notation) 


u^{t + h) — 


We claim that 





(3.10) 


uniformly in N] this would imply p.Op . To prove it we sum up all the estimates coming from Lemma [3.9l 
for different values of k in order to cover the cylinder C. We obtain that there exist C > 0 independent 
of N and t G [ti,t 2 ] such that ||ui(t,/i)^||y < Ch, which implies (|3.10l) . Hence is strongly compact 
in L^iC). Now any compact set in H can be covered by a finite number of open cylinders. To conclude 
we take a countable sequence of compact sets embedded in H whose increasing union exhausts H and 
apply a diagonal procedure. □ 


Corollary 3.11. There exists a subsequenee o/{i2^} whieh eonverges strongly in L^{Qt). 


Proof. We can combine Lemma [3.101 with the uniform bound provided by Lemma [3.5l to use Lebesgue’s 
dominated convergence theorem. Note that the functions are constantly equal to ip outside and 
that Lemma 13.41 applies. □ 
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4 Existence of solutions 

In this section we prove the existence of weak solutions of (j2.ip . 


4.1 Convergence of the approximate solutions 

Lemma 4.1. There are funetions u,u such that the following statements hold (up to extracting a 
subsequence) for iV —oo; 

1) weakly in LP{0 ,T,W^'P{Qo)) n L°°{Qt), 

2) ^ u in L^{Qt) and a.e. in Qt, 

3) —)■ and a.e. in Qt, 

4) ^ u in 

5) u = u + ^ 

Proof. The first statement follows from Lemma 13.51 and Corollary 13.71 The second statement follows 
from Corollary 13.111 To prove the third statement we write 

- V'XQj,\nllLi(Qr) - IIV'lll,°°(QT)IIXQT\n'^ - XQj,\QllLi(Qr) ^ 
as N ^ oo, thanks to Lemma 13.41 It follows that 

U := u -in L^(0). 

Since 11^ —)> 11 by Lemma [TH we get that ^ if a.e. in Qt\(1, so that u is supported on 11. □ 

Recalling Lemma [T5] it follows that, up to a subsequence, —>■ M in LP{Qt) and —)• u in LP{Vl), 

for all 1 < p < oo. 

We now discuss the convergence of the time derivatives. 

Lemma 4.2. There exists A G T>'{Qt) such that, up to extraction of a subsequence, ^ K in 
T^'{Qt)- In fact, A agrees as a distribution over Qt with the time derivative (in distributional sense) 
of the function u defined in Lemma\4(J\ Moreover, given any cylinder C := {ta,tb) x K CC H, there 
holds that l\c € Lp'{ 0,T-,W-^'P'(K)) and u^\c ^ ^\c in Lp'( 0, T; IT-i’P'(iX)). 

Proof. Let us denote by (•, •) the pairing between T>'{Qt) and T>{Qt). Given (f G 'D{Qt), we compute 





u‘^4>t dxdt. 


We may now use Corollary 13. Ill to pass to the limit, so that 



uft dxdt 


lim 

TV—>-00 




up to a subsequence. This shows the first and second statements. 

Our last statement is a consequence of Lemma 13.91 which provides uniform bounds on the time 
derivative over cylinders contained in 11 as in Lemma 13.101 □ 
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Corollary 4.3. Let •= [^i)^ 2 ] x K such that fl 50 = 0. Then u € C{ti,t 2 ', L‘^{K)). 

Proof. Let (f G 'D(Titi,t 2 )- By previous considerations, we know that {4>u)t G L^'{ti,t 2 ',W~^’P'(K)) and 
also {(j)u){t) G ILQ^’^(iL) for a.e. ti < t < t 2 - Using Lemma IHTHl we deduce that (pu G C{ti,t 2 ‘, L'^{K)). 
Being p and K arbitrary, the thesis follows. □ 

Corollary 4.4. There holds that A|^ G V*. 

Proof. Let p G Cf{Ll). Thanks to Lemma 13.91 we have that 

|(A, I?i))v*-v| < ll'?^l|v ('^II^IIlp(o,T,VK1’P(Qo)) ll^llLP(0,r,Lp'(Qo))) ’ 

Our claim follows by a duality argument. □ 

4.2 Recovery of the limit equation 

Our next aim is identifying the limit equation. Let us define 


N-l 

A^{t,x) := Xitk,tk+i)ii)Mtk,x,u'^,Vu'^)xQ(t^){x). 

k=0 

Lemma 4.5. There exists a function A G L^' (Qt)^ such that A‘^ ^ A in L^' (Qt)'^ as N ^ oo, up to 
a subsequence. Moreover A is supported in 0. 

Proof. This follows directly from (j2.3l) and Lemma 13.61 □ 

To identify A we will require a number of auxiliary results. 

Lemma 4.6. Let p be smooth and such that supp p C 0^ 0 0. Given t > 0 we define 

:=- [ {{p{t) - p{s))u{s) ds 

X Jt-T 

(we set /?’’ := 0 when the previous formula does not make sense), being u the function defined in Lemma 
\4.1\ Then p'^ gV for any r > 0 and p'^ ^ Q inV as r —> 0. 

Proof. Since supp p'^ C 0'^ 0 O for small r, we can approximate p'^ in the norm of V by functions in 
C 2 (n) convolving with a mollifying sequence, so that p^ G V. 

Let now K C be an open set such that K C Ll^{t) a.e. t G {ta, h) for some values 0 < ta < tj, < 
T. Thanks to [Ml Ch. 2, Th. 9], we get that p’" —0 in LP{ta,ti,',W^’X{K)) as r —> 0. Covering supp(/) 
with a finite collection of cylinders of the form {ta,th) x K yields the desired result. □ 

Lemma 4.7. Let p be smooth and such that supp p C n 0. Then 

limsup f ( A^-Vu^pdxdt< f ( A-Vupdxdt. (4.1) 

N^oo Jo Jo Jn{t) 


14 


Proof. Let r > 0 and define 


{t) = — [ u{s) ds. 
Jt-T 


u [t] = 
,A 


By multiplying the equation for by {u^ — u^)(j) and integrating by parts we get 

[ iu^ — u^)u^(pdxdt = [ [ {u‘^ — u^) divA‘^(pdxdt 

lo Jn^{t) Jo Jn^(t) 


+ 


fl 


0 Jn^{t) 


lo Jn^(t) 


■ VvPcj) dxdt + 


A^ ■ Vu^cfdxdt — 


Q^{t) 

r 

lo 


A^ -Vcfu^ dxdt 


lo Jn^{t) 


• V(j) vP dxdt :=/ + // + III + IV. 


Let us elaborate on the left hand side of the previous equality. We compute 


/7 

Jo JQ^{t) 


dxdt = 


r [ 

lo Jn^{t) ot 


{u^{t)f 


dxdt 


fl 


{ u ^{ t)y , J , 

--- (ft dxdt —>■ — 


/ / —(ft dxdt 

Jo Jn{t) 2 


as —>■ oo, thanks to Lemma HU Next, we have that 


fl 


0 Jn^{t) 


vAu^fdxdt = 


lo Jn^{t) 

(T 


A "r 
-t- 


u{s) ds dxdt 


t — T 


[ [ — [ (('5^(^) “ dsl dxdt 

Jo I Jt—r J 


fl 

fl 

Jo Jq.‘ 


{<fu)lu^ dxdt — / 

Jo ' 

f{t)u(t) — (f{t — T)u(t — t) 


dxdt 


dxdt — 


lo Jn^{t) 

=: A + B. 

Thanks to our assumptions on (f we have that 

pT p N —1 


lo JQ^(t) 


dxdt 


B = - 


P"^Xit,,t,+t)^tXnit,)dxdt 

lo Ju^it) 


for r small enough. We then pass to the limit in B by Lebesgue’s dominated convergence theorem. 
Indeed, if r is small enough Lemma 14.21 enables to get a.e. convergence of the integrand, domination 
follows as the duality product is uniformly bounded. To deal with the limit of ^ as iV ^ oo we may use 
Lemma I4.1l 4l together with the fact that the incremental ratio is essentially bounded (after Lemma 
[33D . Gathering all the previous and letting —>■ oo, we find that 


fl 


0 


vlut(fdxdt 


f (f{t)u{t) - <f{t - T)u{t - t) 

lo Juit) 


T 


u dxdt — / yo’^A dxdt, 

J Qrp 


15 








which is bounded from below by 


(j){t) - (j){t - t) v?{t) 


Jo Jn{t) 

Letting r —)• 0+ and using Lemma 14.61 we obtain 


/ , 

j Qrp 


dxdt — / K dxdt. 


fl 


v?{t) 

(pt —— dxdt, 


0 Jn{t) ^ 

so that lim inf T-_>.o+ liminfAr_>.oo(/ + // + // + IV) > 0. 

We are now ready to compute the limit of/ + // + /// + IV when N ^ oo. First, we find out that 


Ho¬ 
using Lemmas 14.11^ 41 and 14.51 We also have 

III ^ 


fl 

Jo Jn(t) 


AS/(pu dxdt 


'0 Jn{t) 


AS/vIpdxdt 


as N ^ oo (clearly S/u'^ € IP{Qt)^). Note that S/vI = iS/u)'^ —> S/u in as in the proof of 

Lemma 14.61 Taking limit r —> 0, the above integral converges to 


fl 

Jo Jq 


AS/up dxdt. 


0 JQ.{t) 


Finally, arguing as before we get that 


IV 


'0 JQ.{t) 


AS/pu^ dxdt 


fl 

Jo Jn(t) 


AS/pu dxdt 


as iV —oo, which converges to 

after taking the limit r ^ 0. Hence 

lim sup / / 

n—kx Jo Ju{t) 

and the result follows. 

Lemma 4.8. There holds A{t,x) = A[t,x,u,S/u) a.e. in H. 


A^S/u^pdxdt < 


fl 

Jo Jnit) 


ApS/u dxdt 


□ 


Proof. We use Minty-Browder’s technique. Let 0 < (/> G Cq{Qt) with suppiji C n n, and let 
g G C^{Qt). Thanks to the monotonicity assumption (|2.51) . we have 


jt TV—1 


yi / {A{tk,x,u^,S/u^) — A{tk,x,u^,S/g)){S/u^{t)— S/g)pdxdt >t). 

^0 = \ 
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From Lemma 14.71 we get 


limsup / W / A{tk,x,u^,Vu^)Vu‘^4>dxdt < / / A'Vu4>dxdt. 

N^oo Jo -'o Jn{t) 

We now show that 

/ / A(tk,x,u^,Vg)Vu^4>dxdt ^ / / A{t,x,u,'Vg)Vu4>dxdt, (4.2) 

-'0 Jo Jfl{t) 

as A/" —>■ oo. Indeed, recalling (j2.6p we have 


N-l 

A{t, X, u, Vg) - ^ A(4, a;, Vc/)x[ 4 , 4 +i) 

fc=i 


Af-l 

^ L] ^[h,4+i) (^(1* - tk\) + C\u{t, x)-u^{t, a;)|) \Vg\P-\ 
k=l 


Note that the right-hand side above converges to zero a.e. in O and also in for all p < oo as 

N —>■ oo. On the other hand, Vu‘^ Vu weakly in thanks to Lemma l4.ll which yields (j4.2[) . 

In a similar way we show that 


.tN-i 

/ E / 

Jo 

Finally we obtain that 


A(tk, X, u^, Vg)Vg 4> dxdt 


lo JQ(t) 


A{t, X, u, Vg)'Vg 4> dxdt. 


TV—1 

E 


A(tk,x, u^, Vu^)Vg J) dxdt 


■lO k=\ 

thanks to Lemma 14.51 Summing up, we obtain 

rT 

ln{t) 


lo Jn{t) 


AVg J) dxdt 


f [ {A — A{t,x,u,Vg)){'Vu{t) —'Vg)4>dxdt > 0. 
Jo Jn(t) 


This implies that A = A{t, x, u, Vu) for a.e. (t, x) G suppc/), by means of Minty-Browder’s method (see 
for instance [H Ch. 9.1]). □ 


4.3 Recovery of boundary and initial conditions 

Proposition 4.9. The function u defined in Lemma \4.1\ is a weak solution of vroblem \2.1\ in the sense 
of Definition \2.6l Furthermore, u{t) —>■ uo a.e. as t —>■ 0. 

Proof. Let (f G C^{Qt) with supp</> C fl Q.. We fix a value of A; G {1,..., N" — 1} and test the 
approximating problem in [tfc,t) x Tl{tk) with t < tk+i. That is, 


ln{tu) 

/ u^{tk)4>{tk) dx + / / u^{s)4>sdxds 

-/n(tfc) J^h J0,(tk-) 


u^{t)(j){t) dx + 


• S/(j)dxds 
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for any f € [tk, tk+i)- By adding these contributions from 0 to t l}we get 


f u^{t)4>{t) dx + f f -S/cjidxds 

Jn^(t) Jo Jn^{s) 

uo</'(0) dx + / / u^{s)J)sdxds 

Jo Jn^(s) 

+ V] ( / u^{tk+)(j){tk)dx - [ u^{tk-)(l){tk)dx 


f 

>Q{0) 


(4.3) 


Since supp(/> C we also have 

u^{tk+)4>itk) dx - 


i’{tk)4>{tk) dx - 


u^{tk-)4>(tk)dx 


'n(4)\n(4_i) 


lQ{tk-i)\Q{tk) 


u^{tk-)(t){tk)dx = 0 . 


Thanks to Lemma 14.11 41. converges strongly to u in L^(supp(^). Hence we can pass to the limit in 

(j4.3p and obtain 


ln{t) 


u{t)(p{t) dx + 


0 JQ{s) 


A{t, X, u, Vtt) • S/cpdxds 


ln{o) 


uo4>{0) dx + 



0 Jq,{s) 


u{s)(j)s dxds 


for a.e. 0 < t < T, which holds for any cj) € C^{Qt) with supp(^ C H. This can be stated as 


ut = diYA{t,x,u,Vu) in T>'(H). 

Furthermore, since u € LP{0,T] W^'P{Qq)) and u = 'i/j a.e. Qt\^, we get that u(t) — ip(t) € lTQ^’^(H(t)) 
for almost any t € (0,T). Hence we also recover the boundary conditions at dii} in the limit. 

Let us deal next with the initial condition. Note that for t small enough we have 

/ u{t)(j){t) dx = / uq(I){&) dx + C{(j))t 

Jn{t) Jn(o) 

for some C{4i) > 0. Here we use that we assume condition on the time slicing (and specihcally on 
to = 0) as specihed at the beginning of Section [3l Hence 


lim / u(t)0(t)dx 

Ju(t) 



dx. 


Now let K CC H(0) such that u € (7(0,ti,L^(iL)) for some ti > 0 (which exists as H is Lipschitz). 
Then u{t) converges in Lp‘{K) to some fto as t —)■ 0. This limit uq must agree with the distributional 
limit uq over K. Hence u{t) —)> uq in L^^^(H(0)) as t —>■ 0. In particular we get a.e. convergence to the 
initial condition. Note that this works in the same way for any relatively open subset of d-iQ. 

Finally we justify that u — ip ^V. Once we have shown that the boundary conditions on 9^0 are 
fulfilled, it is easy to construct a sequence rjn belonging to (7^(0) and satisfying \\{u — ip) — r]n\\v —^ 0 as 
n —7> oo. For instance, we may consider G G C^{1R) such that |(7(t)| < \t\, G{t) = 0 if |t| < 1 and G{t) = t 
if \t\ > 2. We also consider pn to be a standard mollifying sequence. Then rj^i = G{npn * {u — ip))/n 
has the desired properties. □ 
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The argument above also shows that, given a cylinder [ti,t 2 ] x K CC fi, the map t i-7> Ujx is 
L^-continuous in [^ 1 ,^ 2 ]- As a consequence, if we fix t > 0 then u(s) —>■ u{t) as s —>■ t a.e. in ^(t). In 
this sense, we can claim that t 1 —>■ u{t) E C(0,T, 


5 Uniqueness of solutions 

We start with a technical result which can been proved as in [281 Proposition 2.6]. 


Proposition 5.1. Let Assumptions [Kl[\2.4\ and E3 be satisfied. Then the following integration by 
parts formula holds: 


rt2 




{ut,v)s + {vt,u)sds = / u{t 2 -)v{t 2 -) dx - / u{ti+)v{ti+)dx, 


(5.1) 


for any 0 < ti < t 2 < T and any n, n E V, where (•, •)* indicates the pairing between W {Ll{t)) and 
Proof of Theorem \2.!A Let ui , U 2 be two solutions of (12. Ij) . Let e > 0 and define 

f • ( \ ( 5|x|^ 2|xp 9\x\^ 4|a;|\ 1 1 , r, 

s,(x) := (-1^ - -y- - 4E + T) 1^1 < g cHm. 

\ sign(x) |x| > 2e 

which is a regularization of the sign function that converges pointwise as e —>■ 0. Note also that we have 
geiui — U 2 ) E LP{0,T, WQ’^{Qfi)). Besides, supp 5 e(ni — U 2 ) lies in the closure of Tl. Then, with a slight 
abuse of notation, ge{ui — U 2 ) = gfiui — U 2 ). 

We pick {4>n\n G T^{Qt) such that 4>n —>■ ge{ui —^ 2 ) strongly in L^’(0, T, ITo^’^(Qo)) and supp (/>„ C Ll. 
Note that the pairing 

((ni — U 2 )t, 4>n)v*-v 

makes sense and is bounded independently of n. Then we substitute (fn in (|2.13l) . On one hand, when 
n —7> 00 we get 

4>niui - U2)t dxdt ge{ui “ U2)iui “ U2)t dxdt. 

J Q. J Q. 

On the other hand, integrating by parts and using (|2.3p . 

(t>n{ui — U 2 )t dxdt = — I ^4>n (A(t, X, ui, Vtti) — Aft, X, U 2 ,Vu 2 )) dxdt 

J D. J D. 

—)■ — — U 2 ) {A{t,x,ui,\/ui) — A{t,x,U 2 ,Vu 2 )) dxdt asn^oo. 

Jn 

Thus, we have shown that 


((ui -U2)t,ge{ui -U2))v*-V = - - U2)V {ui - U2) [A(t, X, Ui , Vui) - Aft, X, U2,VU2)] dxdt. 

Jn 


Using the fact that 


[gfiui - U 2 )]t = g'fiui - U 2 ) ■ (ui - U 2 )t in V'{Tt) 
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and denoting 

( xg'^{x) xG(-2e,2e) 

Peix) := { G C\]R) 

[ 0 1^1 > 2e 

we may argue as before to obtain that 

{[geiui -U2)]t,ui -U 2 )v*-V = - L^\pe{ui - U2)][A{t, X, ui, Vui) - A{t, X, U2,Vu2)] dxdt. 

Jn 

In such a way, 

((ui - U2)t, ge{ui - U2))v*-V + {[geiui - U2)]t,Ui - U2)v*-V 
= - / V(ni - U2)[A{t,x,ui,S/ui) - A{t,x,ui,'Vu2)] {g'eiui — U2) +p[{ui — U2)} dxdt 

Jo. 

— V(ni - U2)[A{t,x,ui,'Vu2) — A{t,x,U2,S/u2)] ^ 2 g'^{ui — U2) + (ui — ^ 2 ) 5 '^ {ui — U 2 )| dxdt 

The first term above is less or equal than zero due to (12.51) and the fact that g'^ +p' > 0, hence we can 
neglect it. As regards the second term, we notice that there is some C > 0 such that 

\9eix)\ < C/e, \xg”{x)\ < C/e Vx G -(2e,2e). 

Then we use (12.6p to write 

2C f 

II < - / X{|«i-« 2 |< 2 e}|V(t 6 l - U 2 )||Vu 2 |^"Vl - U 2 \dxdt 

e JQ 

<AC I I \Viui-U 2 )\\Vu 2 \P~^dxdt := e{e), 


'll 

Jo J\u 


'\u-i—U2\<2e 

which is uniformly bounded with respect to e. In fact this term vanishes in the limit e —)■ 0 given that 
V(tti — U 2 ) = 0 almost everywhere on the set of points such that ui — U 2 = 0. Then, thanks to (j5.ip 
we obtain that 

geiui - U 2 ){T-){ui - U 2 ){T-) dx- [ geiui - U 2 )( 0 )(tti - n 2 )( 0 ) dx < 9{e) 

Jq(o) 


Jn{T-) 

and thus taking the limit e —)■ 0 we find 


/ \ui — U 2 \{T—) dx < 
Jn{T-) 

for any T > 0. This implies our uniqueness result. 


/n(o) 


\ui — 1 ^ 21 ( 0 ) dx 


□ 


Remark 5.2. This proof can be considerably simplified if the operator A does not depend explicitly 
on u, as we can choose ge{x) = x in the previous computations and all the proof boils down to the 
monotonicity property (12.51) . 

Remark 5.3. Let us note that the same uniqueness proof can be extended to the case in which there 
exists a finite number of times to := 0 < < • • • < < t^ '■= T such that ((ti,tj+i) x Qq) n H 

verifies Assumption 12.81 for each i = 0,..., A^ — 1. Namely, the former proof would show that any two 
solutions ui,U 2 with the same initial datum agree on ((0,ti) x Qo) n Taking traces at ti— we find 
that ui = U 2 a.e. on II(ti—). Thus ui = U 2 a.e. on n(ti+) and we can repeat the former uniqueness 
proof to obtain that ui agrees with U 2 on ((ti,t 2 ) x Qo) H ^ and hence on (0,t2) x Qo (~l fl. We can 
continue in this way until we reach uniqueness in the whole of fl. 
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Remark 5.4. We observe that Assumption 12.81 could be replaced by the more general requirement 
that the domain satisfies ()5.ip . In fact, it suffices to have (15.ip with a ”>” instead of and only 
for functions tt,u G V such that uv > 0. 
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